Can OPERA help in constraining neutrino non-standard interactions? 
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We study how much the unique ability of the OPERA experiment to directly detect v T can help 
in probing new, non-standard contact interactions of the third family of neutrinos. We perform a 
combined analysis of future, high-statistics MINOS and OPERA data. For the case of non-standard 
interactions in to transitions we also include the impact of possible DoubleCHOOZ data. In 
all cases we find that the v T sample of OPERA is too small to be statistically significant, even if one 
doubles the nominal exposure of OPERA to 9 x 10 19 pot. OPERA's real benefit for this measurement 
lies in its very high neutrino energy and hence very different L/E compared to MINOS. 



INTRODUCTION 

The confirmation of the neutrino oscillation interpretation of solar and atmospheric neutrino data by 
reactor 1] and accelerator 3] neutrino experiments brings a unique picture of neutrino physics in terms of 
three-neutrino oscillations [4j , leaving little room for other non-standard neutrino properties [5J] . Nevertheless, 
it has long been recognized that any gauge theory of neutrino mass generation inevitably brings in dimension- 
6 non-standard neutrino interaction (NSI) terms. Such sub-weak strength operators arise in the broad class 
of seesaw- type models, due to the non-trivial structure of charged and neutral current weak interactions 0]. 
Similarly, NSI also appear in radiative models of neutrino mass. They can be of two types: flavor-changing 
(FC) and non- universal (NU) and their strength sGf is highly model-dependent but may lie within the 
sensitivities of currently planned experiments. The presence of NSI leads to possibly new resonant effects 
in the propagation of astrophysical neutrinos , 1 1 11 ] and it is interesting to scrutinize their possible 



role in the propagation of laboratory neutrinos. With neutrino oscillation physics entering the precision 
age 



12 



131 ] it becomes an important challenge to investigate the role of NSI in future terrestrial neutrino 
oscillation experiments. 

The interplay of oscillation and neutrino non-standard interactions (NSI) was studied in [141 ] and sub- 
sequently it was shown 15, flfj| that in the presence of NSI it is very difficult to disentangle genuine 



oscillation effects from those coming from NSI. The latter may affect production, propagation and de- 
tection of neutrinos and in general these three effects need not be correlated. It has been shown that 
in this case cancellations can occur which make it impossible to separate oscillation from NSI effects. 
Subsequently it was discovered that the ability to detect v r may be crucial in order to overcome that 
problem [Tjj, though this method requires sufficiently large beam energies to be applicable. Barring the 
occurrence of fine-tuned cancellations, NSI and oscillations have very different L/E dependence. There- 



fore, combining different L/E can be very effective in probing the presence of NSI. The issue of NSI and 
oscillation in neutrino experiments with terrestrial sources has been studied in a large number of publica- 
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32j . In [2|| it was shown that MINOS g|j on 
its own is not able to put new constraints on NSI parameters. On the other hand, in |23j | the combination 
of atmospheric data with MINOS was proven to be effective in probing at least some of the NSI parameters. 
Since matter effects are relatively small in MINOS, its main role in that combination is to constrain the 
vacuum mixing parameters. 

The question we would like to address here is whether the combination of MINOS and OPERA 
can provide useful information on NSI. OPERA has recently seen the first events in the emulsion cloud 
chamber 34j and hence it appears timely to ask this question. The idea is that OPERA will be able to 



detect v T and has a very different L/E than MINOS. Both factors are known to help distinguishing NSI from 
oscillation effects. Clearly, much larger improvements on existing sensitivities are expected from superbeam 
experiments like T2K 35( and NOt/A [361) especially in combination with reactor neutrino experiments like 
or Daya Bay [39( , see Ref . [32j . In this letter we will focus on the simple case where 
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NSI only affects neutrino propagation. 



BASIC SETUP 



Adding NSI into the propagation of neutrinos yields the following evolution Hamiltonian 
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where we have made use of the fact that all e xfl are fairly well constrained and hence are expected not to play 
a significant role at leading order. The effect of e ee is a re-scaling of the matter density and all experiments 
considered here are not expected to be sensitive to matter effects. Hence we will set e ee — 0. Note, that the 
e as defined here, are effective parameters. At the level of the underlying Lagrangian describing the NSI, the 
NSI coupling of the neutrino can be either to electrons, up or down quarks. From a phcnomcnological point 
of view, however, only the (incoherent) sum of all these contributions is relevant. For simplicity, we chose 
to normalize our NSI to the electron abundance. This introduces a relative factor of 3 compared to the case 
where one normalizes either to the up or down quark abundance (assuming an isoscalar composition of the 
Earth), i.e. the NSI coupling to only up or down quark would need to be 3 times as strong to produce the 
same effect in oscillations. Since both conventions can be found in the literature, care is required in making 
quantitative comparisons. 

There are two potential benefits beyond adding statistics from combining the data from MINOS and 
OPERA: First, OPERA can detect v T which, in principle, allows to directly access any effect from e XT . 
Moreover, although the baseline is the same, the beam energies are very different (E) ~ 3GeV for MINOS, 
whereas (E) ~ 17GeV for OPERA. 
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Experiments 



All numerical simulations have been done using the GLoBES software [40j, |4l| . In order to include the 
effects of the NSI we have customized the package by adding a new piece to the Hamiltonian as shown in 
equation[T] We have considered three different experiments: MINOS, OPERA and DoubleCHOOZ, the main 
characteristics of which are summarized in table |TJ 

MINOS is a long baseline neutrino oscillation experiment using the NuMI neutrino beam, at FNAL. It uses 
two magnetized iron calorimeters. One serves as near detector and is located at about 1 km from the target, 
whereas the second, larger one is located at the Soudan Underground Laboratory at a distance of 735 km 
from the source. The near detector is used to measure the neutrino beam spectrum and composition. The 
near/far comparison also miti gate s the effect of cross section uncertainties and various systematical errors. In 



our simulations, based on [42l. |43|. |44| . we have used a running time of 5 years with a statistics corresponding 
to a primary proton beam of 5 x 10 20 per year, giving a total of 2.5 x 10 21 , the maximum reachable value 
reported by the MINOS collaboration. The mean energy of the neutrino beam is (E) ~ 3 GeV. 

The OPERA detector is located at Gran Sasso and gets its beam from CERN (CNGS). OPERA consists 
of two parts: a muon tracker and an emulsion cloud chamber. The latter one is the part which is able 
to discern a u T c harg ed current interaction by identifying the subsequent r-decay. The baseline is 732 km. 
Following 
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42 



45| we assume a 5 year run with a nominal beam intensity of 4.5 x 10 19 pot per year. The 
CNGS neutrino beam has an average energy of (E) ~ 17 GeV . 

Since both MINOS and OPERA have the same baseline we use the same matter density which we take 
constant and equal to its value at the Earth's crust, that is p — 2.7 g/cm 3 . 

Finally, DoubleCHOOZ is a reactor experiment, to be located in the old site of CHOOZ, in France. The 
experiment consists of a pair of nearly identical near and far detectors, each with a fiducial mass of 10.16 1 of 
liquid scintillator. The detectors are located at a distance of 0.2 km and 1.05 km respectively. As considered 
in [46J we assume the thermal power of both reactor cores to be 4.2 GW and a running time of 5 years. The 
neutrinos mean energy is (E) ~ 4 MeV. 

Concerning the neutrino oscillation parameters used to calculate the simulated event rates, we have taken 
the current best fit values given in Ref. unless stated otherwise: 

sin 2 ^ = 0.32, sin 2 6>J, 3 U0 = 0.5, sin 2 0f 3 ue = 0, 

(Am^) truc = +7.6 x 10~ 5 eV 2 , (Am^Y™ = +2.4 x 10" 3 eV 2 , = 0. " 

Note the positive sign assumed for (Am^)'™ which corresponds to the case of normal hierarchy. Since, 
none of the experiments considered here is very sensitive to ordinary matter effects, our results would be 
very similar when choosing as true hierarchy, the inverted one. 



Label 


L 


(E v ) 


power 


trun 


channel 


MINOS 2 (M2) 


735 km 


3 GeV 


5 x 10 20 pot/yr 


5 yr 




OPERA (O) 


732 km 


17 GeV 


4.5 x 10 19 pot/yr 


5 yr 




DoubleCHOOZ (DC) 


0.2 km (near) 
1.05 km (far) 


4 MeV 


8.4 GW 


5 yr 


u e — > v e 



TABLE I: Main parameters of the experiments under study. 
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RESULTS 



Disappearance - Probing NU NSI (e TT ) 



As it has been previously shown in [23j, [24j the presence of NSI, notably e TT , substantially degrades 
the goodness of the determination of the "atmospheric" neutrino oscillation parameters from experiment. 
Indeed as shown in figure Q] our calculation confirms the same effect, showing how the allowed region in the 
sin 2 023-Am^-plane increases in the presence of NSI. 

This figure is the result of a combined fit to simulated OPERA and MINOS data in terms of the "at- 
mospheric" neutrino oscillation parameters, leaving the mixing angle #13 to vary freely. The inner black 
dot-dashed curve corresponds to the result obtained in the pure oscillation case (no NSI). As displayed in 
the figure, allowing for a free nonzero strength for NSI parameters e TT and e er the allowed region grows 
substantially, as seen in the solid, red curve. Intermediate results assuming different upper bounds on |e TT | 
strengths are also indicated in the figure, and given in the legend. One sees that the NSI effect is dramatic 
for large NSI magnitudes. However, such large values are in conflict with atmospheric neutrino data l^ . . 
In contrast, for lower NSI strengths allowed by the atmospheric + MINOS data combination [231 ]. say 
\e Tr \ = 1.5, the NSI effect becomes much smaller. Clearly beam experiments currently can not compete 
with atmospheric neutrino data in constraining e TT . The reason for the good sensitivity of atmospheric data 
to the presence of NSI is the very large range in L/E, especially the very high energy events are crucial in 



constraining NSI 47 1 



In summary, the inclusion of OPERA data helps only for very large values of s TT as can be seen also from 
the first line of tabled]] These large values, however are already excluded by the combination of MINOS and 
atmospheric results [23(. We checked that doubling the OPERA exposure does not change this conclusion. 
The slight improvement by OPERA is exclusively due the sample in the muon tracker and the results do 
not change if we exclude the v T sample from the analysis. The usefulness of the sample stems from the 
very different value of L/E compared to MINOS. These results are not too surprising, since even a very high 
energy neutrino factory will not be able to improve the bound on e TT in comparison to atmospheric neutrino 
data 12611 - 





M2 


O 


M2+0 




90% C.L. 


95% C.L. 


90% C.L. 


95% C.L. 


90% C.L. 


95% C.L. 


Err 


[-10.8,10.8] 


[-11.8,11.8] 


[-10.4,10.4] 


[-11.0,11.0] 


[-8.5,8.5] 


[-9.2,9.2] 




[-1.9,0.9] 


[-2.3,1.0] 


[-2.1,1.4] 


[-2.5,1.6] 


[-1.6,0.9] 


[-2.0,1.0] 


Amli [10" 3 eV 2 ] 


[2.3,4.5] 


[2.2,4.9] 


[2.0,5.0] 


[2.0,5.3] 


[2.3,3.8] 


[2.2,4.0] 


sin 2 623 


[0.08,0.92] 


[0.07,0.93] 


[0.08,0.92] 


[0.07,0.93] 


[0.12,0.88] 


[0.11,0.89] 



TABLE II: 90% and 95% C.L. allowed regions for e TT , e eT , Am^ and sin 2 #23 for different sets of experiments. Each 
row is obtained marginalizing over the remaining parameters in the table, plus #13. The true value for sin 2 29n is 0. 




FIG. 1: Shown is the allowed region in the sin 2 2#23-Am 2 il -plane at 95% CL (2 dof). In this fit #13, £ eT and e TT are 
left free. The different lines correspond to different values for e TT as explained in the legend. 



Appearance - probing FC NSI (e eT ) 

It is well known that, in the presence of NSI, the determination of #13 exhibits a continuous degeneracy [If 
between #13 and e eT which leads to a drastic loss in sensitivity in 6*13. A measurement of only P efl and P^g at 
one L/ E cannot disentangle the two and will only yield a constraint on a combination of #13 and e eT . In this 
context, it has been shown in 17j, that even a very rudimentary ability to measure P^ lT may be sufficient to 
break this degeneracy. Therefore, it seems natural to ask whether OPERA can improve upon the sensitivity 
for Set that can be reached only with MINOS. The latter has been studied in [3] in combination with 
atmospheric neutrinos and on its own in Ref. The result, basically, was that MINOS will not be able 
to break the degeneracy between #13 and e eT and hence a possible #13 bound from MINOS will, in reality, 
be a bound on a combination of e er and #13. 

In table HJ we display our results for a true value of #13 = and no NSI. The allowed range for e eT shrinks 
only very little by the inclusion of OPERA data. As in the case of e TT we explicitly checked that this result 
is not due to the v T sample in OPERA but is entirely due to the different L/E compared to MINOS. Also 
a two- fold increase of the OPERA exposure does not substantially alter the result. 

In order to improve the sensitivity to NSI and to break the degeneracy between #13 and e er it will be 
necessary to get independent information on either e eT or #13. An improvement of direct bounds on e eT is 
in principle possible by using a very high energy v e beam and a close detector, but this would require either 
a neutrino factory or a high 7 beta beam. Both these possibilities are far in the future and will therefore 
not be considered any further in this letter. Thus, we focus on independent information on #13. Reactor 
experiments are very sensitive to #13 but do not feel any influence from e er since the baseline is very short 
and the energy very low which leads to negligible matter effects. This is true for standard MSW-like matter 
effects as well as non-standard matter effects due to NSI We consider here as new reactor experiment 




FIG. 2: Shown are the allowed regions in the sin 2#23-£eT-plane at 95% CL (2 dof). Amfu, #23 and e TT are left free 
in this fit. The solid lines correspond to the combination of MINOS2 and DoubleCHOOZ while the dashed lines also 
include OPERA in the analysis. Each set of lines correspond to different true values for sin 2 2#i3, from left to right: 
0, 0.01, 0.05 and 0.1. 





M2 


O 


M2+0 


M2+0+DC 




90% CL. 


95% CL. 


90% CL. 


95% CL. 


90% CL. 


95% CL. 


90% CL. 


95% CL. 


Ett 


[-10.1,11.0] 


[-11.2,12.0] 


[-10.1,10.3] 


[-10.8,11.0] 


[-7.9,9.0] 


[-8.7,9.6] 


[-5.1,5.3] 


[-5.6,5.8] 




[-4.2,1.3] 


[-4.5,1.5] 


[-4.3,1.5] 


[-5.0,1.8] 


[-3.7,1.2] 


[-4.1,1.4] 


[-0.5,0.4] 


[-0.7,0.5] 


Amii [10~ 3 eV 2 ] 


[2.3,4.6] 


[2.2,5.0] 


[2.0,4.8] 


[2.0,5.2] 


[2.3,4.0] 


[2.2,4.2] 


[2.3,2.8] 


[2.3,2.9] 


sin 2 S23 


[0.09,0.92] 


[0.08,0.93] 


[0.09,0.93] 


[0.08,0.94] 


[0.13,0.90] 


[0.12,0.91] 


[0.24,0.78] 


[0.22,0.80] 



DoubleCHOOZ 



TABLE III: Same as table [TQwith true value sin 2 26>i 3 of 0.1. 
38l | . but for our discussion Daya Bay [39( or RENO [48 : ] would work equally well. In figure 



we show the allowed regions in the sin 29i^-e eT plane for the combinations of MINOS and DoubleCHOOZ 
(red solid curves) and of MINOS, DoubleCHOOZ and OPERA (blue dashed curves) for four different input 
values of sin 2 26> 13 indicated in the plot. As expected, the effect of DoubleCHOOZ in all four cases is to 
constrain the allowed sin 26*13 range. The impact of OPERA, given by the difference between the solid and 
dashed lines, is absent for very small true values of sin 2#i3 and increases with increasing true values. For the 
largest currently permissible values of 6*i 3 ~ 0.16, OPERA can considerably reduce the size of the allowed 
region and help to resolve the degeneracy. In that parameter region a moderate increase in the OPERA 
exposure would make it possible to constrain large negative values of e eT . Again, this effect has nothing to 
do with v T detection and, in this case, is based on the different L/E in i/ e -appearance channel. 
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CONCLUSION 

In this letter we have studied how OPERA can help in improving the sensitivities on neutrino non-standard 
contact interactions of the third family of neutrinos. In our analysis we considered a combined OPERA fit 
together with high statistics MINOS data, in order to obtain restrictions on neutrino oscillation parameters 
in the presence of NSI. Due to its unique ability of detecting v r one would expect that the inclusion of 
OPERA data would provide new improved limits on the universality violating NSI parameter e TT . We 
found, however, that the v T data sample is too small to be of statistical significance. This holds even if we 
double the nominal exposure of OPERA to 9 x 10 19 pot. OPERA also has a sample, which can help 
constraining NSI. Here the effect is due to the very different L/E of OPERA compared to MINOS. This 
makes the OPERA sample more sensitive to NSI. However, the improvement is small and happens in a 
part of the NSI parameter space which is essentially excluded by atmospheric neutrino data. 

We have also studied the possibility of constraining the FC NSI parameter e er . For this purpose it is 
crucial to have a good knowledge of #13. Therefore, we included future DoubleCHOOZ data, since reactor 
neutrino experiments are insensitive to the presence of NSI of the type considered here. Therefore, reactor 
experiments can provide a clean measurement of #13 , which in turn can be used in the analysis of long baseline 
data to probe the NSI. DoubleCHOOZ is only the first new reactor experiment and more accurate ones like 
Daya Bay or Reno will follow. Our result would be qualitatively the same if we would have considered those, 
more precise, experiments, but clearly the numerical values of the obtained bounds would improve. The 
conclusion for e er with respect to the v T sample is the same as before: the sample is very much too small to 
be of any statistical significance. OPERA's different L/E again proves to be its most important feature and 
allows to shrink the allowed region on the sin 2 9i^-e eT plane for large #13 values. Here a modest increase in 
OPERA exposure would allow to completely lift the 9\z-e eT degeneracy and thus to obtain a unique solution. 
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the warm hospitality at IFIC at which parts of this work were performed. This work has been supported by 
Spanish grants FPA2005-01269 (MEC) and ACOMP07/270 (Generalitat Valenciana) and by the European 
Commission RTN Contract MRTN-CT-2004-503369. AEP thanks CERN Theory division for hospitality 
during his stay and MEC for a FPU grant. 
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